Dynamical creation of entanglement by homodyne-mediated feedback 
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For two two-level atoms coupled to a single-mode cavity field that is driven and heavily 
damped, the steady-state can be entangled by shining an un-modulated driving laser on the system 
[S.Schneider, G. J. Milburn Phys. Rev A 65, 042107, 2002]. We present a scheme to significantly 
increase the steady-state entanglement by using homodyne-mediated feedback, in which the driving 
laser is modulated by the homodyne photocurrent derived from the cavity output. Such feedback 
can increase the nonlinear response to both the decoherence process of the two-qubit system and 
the coherent evolution of individual qubits. We present the properties of the entangled states using 
the SO(3) Q function. 

PACS numbers: 42.50.Lc, 42.50.Ct, 03.65.Bz 



I. INTRODUCTION 

The deeper ways that quantum information differs 
from classical information involve the properties, impli- 
cations, and use of quantum entanglement. Entangled 
states are interesting because they exhibit correlations 
that have no classical analog. Any two systems described 
by a pure state that cannot be expressed as a direct prod- 
uct \*f>) — \A)®\B) are entangled (non-separable). For a 
two-qubit system there are four mutually orthogonal Bell 
states, which may be denoted 
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where A and B are two subsystems, |0) and |1) are two 
orthogonal states that could represent the atomic ground 
and excited state, respectively. 

Quantum entanglement is a subject of intensive study 
because it is useful and frequently essential for quan- 
tum teleportation 

mm, 

quantum cryptography [4J, 
uantum dense coding |3| and quantum computation 
Also, entangled atoms can be used to improve fre- 
quency standards Thus, one of the increasing inter- 
ests in this context is to find ways to generate the right 
type of entanglement as well as possible @. There are 
a number of measures, like the concurrence [Tol ]. entan- 
lement of formation |ld| . entanglement of distillation 
ITU, relative entropy of entanglement , and nega- 
tivity [Ulllj], that have been proposed in recent years for 
the purpose of quantifying the amount of entanglement. 

In this paper we address the question: what is the 
maximum amount of steady-state entanglement that can 
be generated in a system consisting of two two-level 
atoms (qubits) collectively damped, with the output be- 
ing measured and fed-back to control the system state? 
This question follows naturally from the study of Schnei- 
der and Milburn |15| which considered the same system 



but without measurement or feedback. In both cases 
the entangled steady state is a mixed state while most 
quantum entanglement algorithms are designed for ideal 
pure states. However, it is usually very hard to create, 
maintain, and manipulate pure entangled states under 
realistic conditions, simply because any system is sub- 
ject to the interactions with its external environment. 
These effects, having their origin in decoherence, may 
turn pure state entanglement into mixed state entangle- 
ment. Therefore it immediately raises the question of 
whether the entanglement can be distilled and used as a 
resource for some quantum communication or computa- 
tion task. 

The paper is organized as follows. The model and the 
corresponding master equation together with the deriva- 
tion of the adiabatic elimination of the cavity mode for 
a two-qubit system are analyzed and the steady-state 
solution is presented in Section [H] The entanglement 
measures and the relationship between purity and entan- 
glement are analyzed in Section IIIII In section IIVI we 
investigate the behaviour of the Q distribution function 
and the density matrix in order to obtain information 
about the entangled states. We discuss our results in the 
concluding section Ivl 



II. MODEL AND MASTER EQUATION 

A quantum computer requires that a set of N two-state 
systems can be prepared in an arbitrary superposition 
state. Each of these systems is said to encode a qubit, 
as distinct from the single bit encoded in a classical two 
state system. 

We consider the case where two qubits are coupled to 
a single cavity mode that is driven and heavily damped. 
The two qubits undergo spontaneous emission into the 
cavity mode, and the fluorescence from the output of 
the cavity is subjected to a homodyne measurement with 
a subsequent feedback to the cavity. Our primary goal 
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of feedback control in this paper is to demonstrate the 
ability of feedback control to increase steady-state en- 
tanglement by counteracting the effects of both spon- 
taneous emission and the measurement back-action on 
the system. The history of feedback control in open 
quantum systems goes back to the 1980s with the work 
of Yamamoto and co-workers |16| , and Shapiro and co- 
workers^^- Their objective was to explain the observa- 
tion of fluctuations in a closed-loop photocurrent. They 
did this using quantum Langevin equations (stochastic 
Heisenberg equations for the system operators) and also 
semiclassical techniques. The latter approach was made 
fully quantum-mechanical by Plimak. For systems with 
linear dynamics, all of these approaches, and the quan- 
tum trajectory approach of R.efs . [TsL l2(Tl | . are equally easy 
to use to find analytical solutions. The advantage of 
Wiseman and Milburn's (quantum trajectory) approach 
to quantum control via feedback is for systems with non- 
linear dynamics, as we will discuss in this paper. 

To obtain a master equation for the two-qubit system, 
we follow the analysis of Schneider and Milburn's paper 
[l5j for the steady-state of a system of two qubits in- 
teracting simultaneously with a driving laser. They have 
shown that the steady state is entangled. In this work, 
we expand their analysis to include feedback modulation 
of the amplitude of the driving laser. 

A schematic diagram of the apparatus is shown in 
FigH To ensure that the two atoms see the same phase 
and amplitude of the cavity mode, we need to locate them 
at points of the standing wave in which the two atoms 
are separated by an integer number of wavelengths. The 
dipole-dipole interaction between the atoms can thus be 
ignored in the following calculation. We assume a strong 
coupling g between the atoms and the cavity mode. The 
fact that the coupling is the same for both atoms means 
that they are indistinguishable. This leads to interference 
in their damping via the cavity, as we shall see. Such 
a proposed approach can be realized experimentally in 
a cavity QED system. In this paper we reproduce the 
results in |15| and show that one can increase the steady- 
state entanglement by using feedback-modulation of the 
laser that drives the cavity mode. 

It is assumed that that the laser interacts with the 
two atoms simultaneously, forcing each atom to undergo 
Rabi oscillations at the same frequency and phase. We 
can thus define the total angular momentum operators 
for the two qubits: 

J ± = &f + *±. (2) 
= J x ±iJ y , (3) 

where 

a = a x ± lay (4) 

In terms of these operators, the cavity-atom interaction 
Hamiltonian is: 

H= g -[J-b+ + j + b] (5) 




FIG. 1: Diagram of the experimental apparatus. The laser 
beam is split to produce both the local oscillator /3 and the 
field «o which is modulated using the homodyne current I(t) 
derived from the damped cavity. 

where here the annihilation operator b + describes the 
cavity mode and g is the coupling constant between the 
qubits and the cavity. Since this is symmetric in the 
atoms, it is natural to use the angular momentum states 
to describe thetwo two-level atoms. The most interesting 
dynamics occurs in the j = 1 subspace. In terms of the 
individual atomic levels, the three states for j = 1 are 

|l) = |e) 1 |e) 2 , \2) = -j=(\g)i\e) 2 +\eh\9h), |3> = l<?)i|<?>2. 

(6) 

The fourth state, i.e. the j = subspace, is 

|* i =o) = ^(|»)i|e> a -|e) 1 | fl ) 2 ). (7) 

This latter ssubspace will not change under any of the 
transformations we perform, as will be seen below when 
we analyze the adiabatic elimination of the cavity mode. 

A. Homodyne Detection 

For simplicity in explanation of homodyne detection, 
let us now consider a system with a single two-level atom 
and subject the atom to homodyne detection. We as- 
sume that all of the fluorescence of the atom is collected 
and turned into a beam. Ignoring the vacuum fluctua- 
tions in the field, the annihilation operator for this beam 
is y/jcr, normalized so that the mean intensity ^(a^a) 
is equal to the number of photons per unit time in the 
beam. This beam then enters one port of a 50:50 beam 
splitter, while a strong local oscillator (3 enters the other. 
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To ensure that this local oscillator has a fixed phase rela- 
tionship with the driving laser used in the measurement, 
it would be natural to utilize the same coherent light field 
source both as the driving laser and as the local oscilla- 
tor in the homodyne detection. This homodyne detection 
arrangement is as shown in Fig. ^ 

Again ignoring vacuum fluctuations, the two field op- 
erators for the light exiting the beam splitter, b\ and b 2 , 
are 

bk= [V7^-(-l) fe /?]/V2- (8) 

When these two fields are detected, the two photocur- 
rents produced have means 

h = {W - (-l) fe (V7/3^ + yfytrp) + l^cx) /2. (9) 

The middle two terms represent the interference between 
the system and the local oscillator. 

Equation @ gives only the mean photocurrent. In 
an individual run of the experiment for a system, what 
is recorded is not the mean photocurrent, but the in- 
stantaneous photocurrent. This photocurrent will vary 
stochastically from one run to the next, because of the ir- 
reducible randomness in the quantum measurement pro- 
cess. This randomness is not just noise, however. It is 
correlated with the evolution of the system and thus tells 
the experimenter something about the state of the sys- 
tem. The stochastic evolution of the state of the system 
conditioned by the measurement record is called a "quan- 
tum trajectory" [2(j. Of course, the master equation is 
still obeyed on average, so the set of possible quantum 
trajectories is called an unravelling of the master equa- 
tion j2jj. It is the conditioning of the system state by 
the photocurrent record that allows feedback control the 
system state at the quantum level. 

The ideal limit of homodyne detection is when the lo- 
cal oscillator amplitude goes to infinity, which in practi- 
cal terms means \f3\ 2 ^S> 7. In this limit, the rates of the 
photodetections go to infinity, and thus each photodetec- 
tor produces a continuous photocurrent with white noise. 
For our purposes the only relevant quantity, suitably nor- 
malized, is the difference between the two photocurrents 

Ell 

= m = + 6 ^)c (*) + 

(10) 

A number of aspects of Eq. I|1U|I need to be explained. 
First, $ = arg/3, the phase of the local oscillator (de- 
fined relative to the driving field). Here we set $ = 0. 
Of course, all that really matters here is the relation- 
ship between the driving phase and the local oscillator 
phase, not the absolute phase of either. Second, the sub- 
script c means conditioned and refers to the fact that if 
one is making a homodyne measurement then this yields 
information about the system. Hence, any system av- 
erages will be conditioned on the previous photocurrent 
record. Third, the final term £(i) represents Gaussian 



white noise, so that 

£{t)dt = dW(t), (11) 

an infinitesimal Wiener increment defined by pl| 

[dW{t)f = dt, (12) 
E[dW(t)] = 0. (13) 

B. Adiabatic elimination of the cavity mode 

Let us now come back to the two-qubit system. The 
cavity mode is uninteresting as for the high levels of 
damping its behaviour is slaved to the driving, so it will 
then be adiabatically eliminated [lH I22I |23| in our first 
calculation, resulting in a master equation followed by 
the density operator p, where p only includes the two 
qubits. 

The complete master equation describing the system 
pumped by an unmodulated driving laser, including the 
cavity mode, is described by a density operator lo given 
as follows. 

Cj = --a[b-b\u]-i^[J-b+ + J + b,oj] 

+71D [eri] u + 7 2 X> [cr 2 ] lu + j p T> [b] lu. (14) 

Here the annihilation operators b describe the cavity 
mode. 71 , 72 and 7 P are the coefficients of damping for 
the two qubits and b mode respectively. I? is a super- 
operator defined as V = T>[A)B = ABA'' - {A^A,B}/2 
for irreversible evolution. The cavity mode is heavily 
pumped and damped. 

The adiabatic elimination of cavity mode is done first 
by displacing the density operator u, and the master 
equation describing its evolution, to zero in the b mode. 
We assume that 7 P is sufficiently large that the state stays 
close to an equilibrium coherent state with amplitude 

(3 = a/ 7p . (15) 

The displacement operator 

D b = e f3(b+ -V (16) 

is used to carry out this transformation. The new density 
operator is v = Db{—(3)coDb(f3). Applying this operation 
to the original master equation in Eq. (|14fl gives the 
new master equation for v in which the b mode is of zero 
average amplitude. This is 

v = Cv - i|[(J + 6 + J~b + ), v\ + lp V [b] v (17) 

in which all the terms involving only the qubits in the 
superoperator C, defined as 

Cu = ~ii^[{J + + J-), v] + 71 2? M v + 12V [a 2 ] v(18) 
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Since the amplitude of mode b is small, a partial ex- 
pansion of the density matrix v in terms of the b mode 
number states need only be carried out to small photon 
numbers. Therefore 

y = Po|0)(0| + (/7 1 |l)(0|+H.c.)+p 2 |l>(l| 

+ ( P2 '|2)(0|+H.c.) + O(A 3 ), (19) 

where A = g/"/ P is a very small number and |0) is the 
cavity vacuum state. This is substituted into the master 
equation Eq. ijTTI) which is expanded and terms multiply- 
ing equal sets of b mode number projectors are gathered 
together to get a set of four equations. Terms of greater 
then second order are neglected. These equations are: 



Po 


= £po - 


■9\r+ 

Pi 


-p\J }+l P P2 




Pi 


= Cpi - 


i^V Po 


+ V2J+ V - j- V] - 


7p 
- Y Pi 


Pi 


= £p2 - 


£\j + pi 


- PiJ +2 } - ■y p p 2 




pi' 


= Cpi- 


-i 9 -[V2J- 


>i] ~ IpPi 


(20) 



Now we make the assumption that both p\ = and p 2 ' — 
so that by using the second and fourth of Eq. I|20|) , the 
values of p\ and p' 2 are found to be 



Pi = — [J Po - J Pi] 

Ip 



Pi 



V2l P 



J Pi, 



(21) 



where pi = 0(— ) and p 2 — 0( — ). These are then 
substituted into the first and third equation of set Eq. I|20l) 
which become 

g i 

po = Cp Q -- — [J + Jp + poJ + J -J + J P2 

-p 2 J + J~] +1 V P2, 

a 2 

pi = C P 2 + — [J- P0 J + - J~p 2 J + ] - 1 P P2- (22) 
7p 



Adding these two equations together and noting that 

2 

p 2 = 0(— ), then neglecting the p 2 terms gives for the 
final master equation of the system. 



p = C P +^-V [J~]p, 

Ip 



(23) 



The coefficient 3— describes the strength of collective 
damping of the two-qubit and the cavity b mode and 
will be called 7. 

The term Cp is expanded to give the master equation 

as 

P = -ii^l(J + +J-),p}+JiV p+ l2 V [a 2 ] p+-/V [J-] p 

(24) 



This is the super-fluorescence master equation [24j. The 
approximation that the timescale imposed by collective 
decay rate 7 P is greater than the timescale of the two 
qubits evolution imposed by single decay rate 71 , 72 , that 
is 7i,72*C7, results in the following master equation 



■ 9 a \i T+ 

p= - l w P [{J 



J-),p]+ 7 V [j-] 



p- 



(25) 



Following the method of Ref . jig , the stochastic master 
equation (SME) conditioned on homodyne measurement 
of the output of cavity is 



dp c = -i[H,p c (t)]+dt^V[J-]p c (t) 
+^dW{t)H[J-] Pc . 



(26) 



The homodyne photocurrent, normalized so that the de- 
terministic part does not depend on the efficiencyis 



C. Dynamics with Feedback 



(27) 



If we now add dynamics with feedback from feedback 
Hamiltonian H fh = I(t)F, where I(t) = (J+ + J~)(t) + 

£(t)/VV, we can S et the SME 

dp c = dt^f'D[J~]p c — idt[H a , p c ]—idt[F, —iJ~p c + ip c J + ] 

+dt-V[F]p c +dW(t)n[-iJ 7 J~ - iXJJpc- (28) 

7 

Here H a — aJ x , F — \J Xl J x = J~ + J + and a, A are 
driving and feedback amplitude respectively. This cor- 
responds to having a feedback-modulated driving laser. 
This is also an Ito stochastic equation, which means that 
the ensemble average master equation can be found sim- 
ply by dropping the stochastic terms. 

Therefore from the above master equation in j = 1 
subspace, we can write down the equation of motion for 
the components of the 3x3 density matrix of the state of 
the system, taking into account that Tr(p) — 1 and that 
p is Hermitian. In order to get the steady state solution 
of the above master equation Ea. (|28|l . we define 

Xij = (i\J x \j) + (j\Jx\i), 



Vij = (i\Jy\j) + (j\Jy\i), 
Zij = (AJz\j) + (j\J z \i)- 



(29) 



Then the differential equations for xtj , yij , Zij are found 
to be 

A 2 

±12 = -27x12 -2AH x 23 + V2ayi 3l 

7 

A 2 

±13 = (-7 2\)x 13 - V2ayi2 + V2ay 23 , 

7 

A 2 

±23 = (27 H )xu ~ {l 2 + A 2 )x 23 + \/2ayi 3 , 

7 



5 



5A 2 

2/12 = —a\l2x\z + (27 + 6A H )yi 2 - 2v2ayi 3 , 

7 



2A 2 



2/13 



-V2axi2 + V2ax23 + (7 + 2A H )j/i 3 , 



3A 2 



J/23 = V2aa;i3 + (-27 6A)yi2 

7 

5A 2 

+ (7 + 4A H )t/23 + 2V2a2i2 - 2V2azi 3 , 

7 



212 



213 



(-2A- 
v 3 



3A 2 , /— /- 87 

)xi3 - 2V2a?/i2 + av 2y 23 irZvi 

7 3 



4A 2A 2 
3 7 



^13 



2 7 4A 
~3 3~ : 



/- r 27 4A 

2Axi 3 - v2ayi2 - v2aj/ 23 + (-g- + "y) 2 ^ 

,47 8A 2A 2 4 7 8A 

3 3 7 3 3 



(30) 



Here we have ignored z 23 since we require only 8 param- 
eters. The steady-state solutions are 





A 




B 


C 




£> 


X12 = 


~s> 


X13 = 


S' 


X23 = S 7 ' 


2/12 = 






E 




F 


G 






yi3 = 


5' 


J/23 = 


5' 


Zl2 = 5' 


^13 = 


1?' 



(31) 



Here 



S = 24a 4 7 4 + 2 7 8 + 267 7 A + 1437 6 A 2 + 432 7 5 A 3 

+789 7 4 A 4 + 926 7 3 A 5 + 726 7 2 A 6 + 352 7 A 7 + 96A 8 
+2a 2 7 2 (4 7 4 + 22 7 3 A + 65 7 2 A 2 ) + II67A 3 + 60A 4 , 

(32) 

and 

A = 

B = 2 7 (7 + 2A)(4a 2 7 2 (7 2 +3 7 A + A 2 ) + 

A 2 (2 7 4 + 14 7 3 A + 33 7 2 A 2 + 32 7 A 3 + 16A 4 )) 
C = 

D = 2\72q7 2 (7 + 2A)(4q! 2 7 2 + A 2 (57 2 + 207A+16A 2 )) 
E = 

F = 2\/2Q7 2 (7 + 2A)(4a 2 7 2 + 27 4 + 147 3 A + 

377 2 A 2 + 447 A 3 + 447 A 3 + 16A 4 ) 
G = 7( 7 + 2A)(4a 2 7 2 ( 7 2 +37A + A 2 ) + 

A 2 (2 7 4 + 14 7 3 A + 33 7 2 A 2 + 32 7 A 3 + 16A 4 )) 
H = 7(7 + 2A) 3 (4a 2 A 2 + 2 7 4 + 

147 3 A + 337 2 A 2 +327A 3 + 16A 4 ). (33) 

III. ENTANGLEMENT AND PURITY 

Now that we have the steady-state solution of the mas- 
ter equation, we can determine if feedback-modulated 
driving can increase the steady-state entanglement when 



Concurrence for different 
driving and feedback (a) 



Concurrence when feedback is off (b) 




FIG. 2: Plot of concurrence and purity vs driving and feed- 
back amplitude. 



the collective decay is considered. Let us now specify the 
measures which we will be using to characterize the de- 
gree of entanglement of a state. As we mentioned before, 
there are several measures of entanglement. In this pa- 
per since we have two qubits system and we choose the 
concurrence [T(tI | as a measure for it. For a mixed state 
represented by the density matrix p, the " spin- flipped" 
density operator, which was introduced by Wootters ^(| > 
is given by: 



P = {(Ty 8> (Ty) P (a y ® (Ty) . 



(34) 



where the bar of p denotes complex conjugate of p in the 
basis of {\gg), \ge), \eg), |ee)}, and a y is the usual Pauli 
matrix given by 



-i 

1 



(35) 



In order to work out the concurrence, we need to de- 
termine the square root of the eigenvalues Ai,A 2 ,A 3 ,A4 
of the matrix pp and sort them in decreasing order, i.e., 
Ai>A 2 >A 3 >A4. It can be shown that all these eigenval- 
ues are real and non-negative. The concurrence C of the 
density matrix p is defined as 



C(P) 



A 4 ,0 



(36) 



The range of concurrence C is to 1 . When C is nonzero 
the state is entangled. The maximum entanglement is 
when C = 1. 

Fig. is the plot of concurrence and purity vs driv- 
ing and feedback amplitude. As shown in Fig [21 a2, we 
can get certain amount of entanglement in the steady- 
state of a unmodulated driving system ^J. In this case 
the maximum concurrence is about 0.11, with appropri- 
ate choice of the driving amplitude a = 0.387~ 1 . The 
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steady state is independent of the initial state as long 
as it is in the symmetric subspace, such as \g)(g\- The 
coherent evolution alone is not able to produce any en- 
tanglement for an initially unentangled state, as it only 
consists of single qubit rotation with no coupling between 
two qubits. Therefore the steady-state entanglement is 
due to the common cooperative decoherence coupling to 
the cavity environment acting together with the coherent 
evolution [Tlj . 

After including feedback onto the amplitude of the 
driving on the atom, proportional to the homodyne pho- 
tocurrent, we see that feedback is remarkable as the 
steady-state concurrence has been improved from 0.11 
to 0.31 as shown in Fig[2al, with appropriate choice of 
driving amplitude a — ±0.47 _1 and feedback amplitude 
[3 = -O.87- 1 . 

The gain of the steady state entanglement comes at the 
price of a loss of purity, as shown in Fig. [21 (b). There 
are a number of measures that can be used for the degree 
of purity, for example, the von Neumann entropy given 
by S = — Tr[plnp], and the trace sqared of the density 
matrix. In this paper we choose the measure of purity 
given by 

r 2 = \{^[P 2 \ ~ \)- (37) 
From the above equations, we find that 

r = g + g [ X 12 + X 13 + X 23 + 2/l2 + 2/l3 + X 23 + Z 12 

+ i(22 13 -z 12 ) 2 ]. (38) 

The minimum purity in 3 x 3 subspace is obviously g. 
Note that Eq(36) measure is linear in Tr[p 2 ] with min- 
imum and maximum 1. The minimum of is only 
attainable in 4 x 4 Hilbert space. 

To gain further insight into the purity of the steady- 
state when it is entangled, we may look at the position of 
the steady-state located in the concurrence-purity plane. 
We begin by choosing a series of driving and feedback am- 
plitudes, and determine their corresponding purity and 
entanglement. In Fig[3| we display these results with 
dots, representing the case of feedback-modulated driv- 
ing, and stars, representing the case of unmodulated driv- 
ing. Obviously, the feedback mechanism leads to a no- 
ticeable increase of entanglement, though it results in a 
less pure state. The continuous curve in Fig|3] represents 
the maximally- entangled mixed states, states with the 
maximal amount of entanglement for a given degree of 
purity, or in other words, states with the minimum pu- 
rity for a given concurrence. The concurrence and the 
purity of a maximally entangled mixed state satisfy the 
following equation pj 

r 2 = l^[4 3 (2-3. 9 )-C7 2 ], (39) 

where 

_ f C/2 for C > 2/3 
9 ~ \ 1/3 for C < 2/3 




0.2 0.4 0.6 0.8 1 
purity 



FIG. 3: Plot of concurrence vs purity. The dots represent- 
ing the case of feedback-modulated driving laser with two 
adjustable parameters a and A, and the stars the case of un- 
modulated driving only with one adjustable parameter a. The 
continuous diagonal curve tands for the maximally-entangled 
mixed states, states with the maximal amount of entangle- 
ment for a given degree of purity 

In reference Q, the degree of the mixture of a state is 
defined by linear entropy Sl- The purity r 2 defined in 
this paper is related to Sl by r 2 = 1 — Sl- 

IV. Q FUNCTION AND DENSITY MATRIX 

To gain further knowledge about the nature of the 
steady-state when it is maximally entangled, we look at 
the density matrix and Q function. Since the density ma- 
trix elements are complex numbers, we plot in Fig^Jthe 
modulus of the matrix elements. We use the unentangled 
basis states 

IVh) = |g)ils>2, (40) 

|^2> = |5>i|e>2, 
1^3) = |e)i|g>2, 
|V^4> = |e)i|e) 2 . 

because a separable basis is needed to discuss entangle- 
ment. 

To define the Q function we need atomic coherent state 

10 ism 

10 = (l + l(C)| 2 )- J exp[CJ+]|i,-i>. (41) 

The atomic coherent state is the closest quantum me- 
chanical states to a classical description of a spin system. 

The Q function is a positive distribution function, 
which is defined as 

Q(M) = (j.MIp&M) 
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(j + m)!(j - m)\ 



E 

m=—j 



cos j+m (-)sm 3 - m (-)e-" n ' 1 



\h m) , 
(42) 



where the state \j, 8, (ft) is the spherical representation of 
the coherent state |C) where 6 and </> are the standard 
spherical polar coordinates defined by £ = arctan^e -11 ^. 
The magnitude of the Q function in a particular direction 
(6 and 4>) is represented in FigQJby the distance measured 
from the origin. Since the Q function is the projection of 
the density matrix into a coherent state , the Q function 
does not give more information than the density matrix. 
However, the advantage of the Q function is that it gives 
a more intuitive view of where the state is located. 

We plot the steady-state Q function and density matrix 
in Fig0] There are a number of points that need expla- 
nation. First, when there is only the unmodulated driv- 
ing laser shining on the system, the steady state having 
the most entanglement is mainly confined to the ground 
state Q function and the upper state is almost unpopu- 
lated. However, when the driving amplitude is modified 
by feedback, the upper state becomes very well popu- 
lated and the entanglement is greatly increased. This 
is not surprising, as more population in the upper state 
enhances nonlinearity both in the decoherence coupling 
and *coherent* evolution of the two-qubit system. 

Second, in the particular separable 4x4 basis , the 
off-diagonal elements of the density matrix represent co- 
herence and the presence of coherence is a necessary con- 
dition for the creation of entanglement. To illustrate this, 
we plot three Bell states which are the maximally entan- 
gled states in FigQ] We see that the off-diagonal terms 
are present in these maximally entangled states. In con- 
trast, we also plot three non-entangled states, the iden- 
tity, the ground state and the excited state. We see that 
none of the off-diagonal elements appear. 

Third, when there is only the unmodulated driv- 
ing laser shining on the system, the density matrix of 
the most entangled steady-state looks most similar to 
that of superposition of \gg) and some of a Bell state 
^j(|e<?) + ls e ))i while when the feedback modulation 
is switched on, the steady-state density matrix looks 
more like a mixture of a non-maximally entangled state 
a\eg) + b\ge) and |ee) + \gg). 



V. DISCUSSION 

To summarize, entanglement between two qubits can 
be created dynamically by driving and coupling them 
to a heavily damped cavity mode. When there is only 
the unmodulated driving laser shining on the system, the 
maximum steady-state concurrence (a measure of entan- 
glement) is 0.11 |l5j. In this paper we have re-derived 
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FIG. 4: Plots of the Q function and its projections on each 
of three planes, and the absolute values of matrix elements of 
the density matrix are shown for each of eight different states, 
as indicated in the figure. The magnitude of the Q function in 
a particular direction is represented by the distance measured 
from the center of the Q function. The diagonal elements of 
the density matrix are |ee)(ee|, \eg)(eg\, \ge){ge\, \gg){gg\ from 
left to right respectively. 



these results and constructed a scheme to increase the 
steady state entanglement by using homodyne-mediated 
feedback, in which the driving laser is modulated by the 
homodyne photocurrent derived from the cavity output. 

An analytical form for the steady-state solution of 
the master equation with feedback was derived and was 
used to show that the amount of the maximum steady- 
state concurrence has been increased from 0.11 to 0.31. 
The properties of the entangled state were also studied 
through the discussion of the Q function and the density 
matrix. The important point here is that the feedback 
scheme can lift the steady-state from the ground state 
towards the excited state. Indeed with such a feedback 
scheme the most entangled state much closer to the Bell 
states, which are maximally entangled states. The details 
about how the feedback mechanism dynamically changes 
the position of the steady-state when it is maximally en- 
tangled are still largely unknown and require further in- 
vestigation. Open questions such as "can the most en- 
tangled steady-state be realized experimentally?" , "what 
happens if individual decay of each atom cannot be ig- 
nored in the calculation" , "how does one increase both 
entanglement and the purity of a mixed state" , "to what 
extent can the system be treated semi-classically" and 
"what happens for a higher number of atoms" are sub- 
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jects for future exploration. 
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